Longitudinal vortices in granular flows 
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We present a new instability observed in rapid granular flows down rough inclined planes. For 
high inclinations and flow rates, the free surface of the flow experiences a regular deformation 
in the transverse direction. Measurements of the surface velocities imply that this instability is 
associated with the formation of longitudinal vortices in the granular flow. From the experimental 
observations, we propose a mechanism for the longitudinal vortex formation based on the concept 
of granular temperature. 
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In fluid mechanics, flow instabilities lead to the for- 
mation of patterns which strongly affect the dynamics, 
yielding coherent structures and controlling the transi- 
tion to turbulence p|. For granular flows, the relevance 
of these hydrodynamic concepts and the existence of a 
similar scenario are open issues. Although granular me- 
dia can exhibit fluid- like properties, they do not behave 
like classical fluids [3] . The dissipative nature of grain 
interactions and the overlap between grain scale and flow 
scale are fundamental differences with classical fluid flows 
which lead to difficulties in our research for a hydrody- 
namic description of granular flows [Q . One way to better 
understand the specificities of granular flows is to study 
the instabilities that can develop in these media. Pattern 
formation has been investigated in vertically vibrated 
granular layers where the non trivial balance between the 
injected energy and the dissipation due to collision leads 
to a rich phenomenology H || . The instability of simple 
shear flows leading to the formation of clusters has also 
been studied numerically and theoretically @, . In con- 
trast, few studies concern the stability of granular flows 
along slopes which are the paradigm for geophysical situ- 
ations. The formation of roll waves has been reported || 
and a fingering instability observed during the propaga- 
tion of a granular front has been described . In this 
letter we present a new instability that spontaneously 
generates longitudinal vortices (i.e. with axes parallel to 
flow) in rapid granular flows down rough inclined planes. 
Although such structures are common in fluid mechanics 
(Gortler vortices |Ll], [l^, streaks in turbulent boundary 
layers |l3| , they have never been observed in a granu- 
lar flow. From the experimental observations, we propose 
a mechanism for the longitudinal vortex formation which 
is specific to rapid granular flow and relies on the concept 
of granular temperature Jl5| . 

The experimental set-up consists of a rough inclined 
plane with a reservoir containing the granular material 
(Fig. 1). The plane is a glass plate (1.3 m long, 0.3 
m wide) made rough by randomly gluing one layer of 
particles onto its surface. The glued particles are the 
same as the ones used for the flow and cover about 70 
% of the plane. The side walls are glass plates and have 
no influence on the flow considering the large width of 
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FIG. 1: Sketch of the experimental set-up. The three pictures 
correspond to top views of the free surface lit from the side 
at three different locations along the slope. 9 — 41°, h g =13 
mm. The sheet laser light is used to measure the surface 
deformation (see picture in Fig. 2). 



the channel. The results presented here have been ob- 
tained using sand with a mean diameter d =0.25 ± 0.03 
mm. Experiments carried out with 0.8 mm sand and 
0.5 mm monodispersed glass beads give the same results. 
The amount of particles in the reservoir is always large 
enough to ensure a constant flow rate during the whole 
experiment. Therefore, the only two control parame- 
ters of the experiment are the inclination of the plane 
9 and the opening of the gate h g . This configuration has 
been widely used to study the rheology of granular flows 
|l8[ ], instabilities |l(| and avalanches dynamics 
and provides a simple laboratory model for geologi- 
cal granular flows down slopes. 

Here we study a new instability observed when the in- 
clination 9 and the opening h g are large (Fig. 2). In 
this regime, the granular material flowing out from the 
reservoir accelerates along the slope while the thickness 
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of the granular layer decreases. At a certain distance 
from the outlet (from 0.4 m to 1.3 m depending on the 
flow conditions), a regular pattern develops and longitu- 
dinal streaks parallel to the flow direction are observed 
(Fig. 1). In the pictures of Fig. 1, the plane is lit from 
the side. The bright and dark stripes are thus the signa- 
ture of a periodic free surface deformation. A sheet laser 
light projected at a low incident angle provides a pre- 
cise measurement of the local thickness (Fig. 2). About 
0.2 m downstream of the initiation of the instability, 
the pattern reaches a "saturated" state: the deforma- 
tion amplitude Ah between crests and troughs and the 
pattern wavelength A remain constant along the slope. 
In this saturated state the mean flow also becomes uni- 
form: mean thickness h, mean velocity and mean volume 
fraction no longer vary along the slope. The wavelength 
A seems to scale with the mean layer thickness: A ~ 3 h. 
However, it is difficult to significantly vary h. In the sat- 
urated state, h only varies between 2.1 and 2.6 mm when 
varying the height of the gate h g or the inclination 6 over 
the range indicated in Fig. 2. 
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FIG. 2: Phase diagram in (8, h g ) plane: the instability is 
observed in the hatched zone. The picture is a close-up of 
the free surface in the saturated regime (8 — 41° and h g =13 
mm). The corrugated laser line provides the thickness pro- 
file measurement: mean thickness h ~2.5 mm, deformation 
amplitude Ah ~0.5 mm, wavelength A <~7.5 mm. 

The streaked pattern described above is not station- 
ary but slowly drifts in the transverse y-direction with 
a phase velocity small compared to the chute x-velocity 
(the drift is about 1 cm/s while the chute velocity is of 
the order of 1 m/s). This drift leads to a complex non- 
linear spatio-temporal evolution including annihilations 
and creations of streaks, which is beyond the scope of this 
study. In this letter we focus on the flow structure and 
try to understand the underlying instability mechanism. 

We have first investigated the grain motions once the 
pattern is fully developed, by measuring the velocity field 
at the free surface. For this purpose, black tracers of the 
same size of the bulk particles arc mixed in the material. 
The flow being rapid, a fast CCD (charge-coupled device) 
camera (4500 frames/s, Kodak HS4540 ektapro) is neces- 
sary to capture the particles displacement at the surface. 



The movies of particles displacement are analyzed us- 
ing a Particle Image Velocimetry method |2C}| . Typical 
longitudinal and transverse velocity profiles V x (y) and 
V y (y) are given in Fig. 3a. The instability induces spa- 
tial velocity modulations which are correlated with the 
surface deformation. First, one observes that the longi- 
tudinal particle velocity V x is no longer uniform across 
the bed but becomes greater in the troughs than in the 
crests with the modulation of AV X being about 20 cm/s. 
The second result is that the instability also induces pe- 
riodic modulations of the transverse velocity V y (AV y ^4 
cm/s): particles no longer follow the bed slope but also 
experience lateral motions. The important observation 
is the 7r/2 phase shift between V x and V y : in the troughs 
V x is maximum whereas V y vanishes. This implies that 
particles move from the crests towards the troughs as 
sketched in Fig. 3b. 
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FIG. 3: a, Longitudinal and transverse free surface velocity 
profiles V x (y) and V y {y) {9 = 41°, h 9 = 13 mm). The verti- 
cal solid (respectively dashed) line indicates the position of 
troughs (respectively crests). To get these measurements, a 
30 mmx30 mm area of the free surface (at x = 1.1 m) is im- 
aged with a 4500frames/s fast CCD camera. The PIV method 
gives the velocity field V x (x,y) and V y (x,y) from two time- 
successive pictures with a spatial resolution of 0.3 mm. V x (y) 
and V y (y) are obtained by averaging the velocity field along 
the x-axis and by averaging over 50 successive pictures (11 
ms). The error bars represent the dispersion of the measure- 
ments, b, Sketch of the particule trajectories showing the 
longitudinal vortices. 

This two-dimensional picture of the surface flow im- 
plies a three-dimensional particle motion and the pres- 
ence of longitudinal vortices in the bulk. The mean flow 
being uniform, mass conservation indeed requires that 
particles move upwards at the crests and downwards at 
the troughs. The longitudinal vortices are then counter- 
rotating with one wavelength A of the wavy surface cor- 
responding to a pair of vortices as sketched in Fig. 3b. 
This 3D structure is reminiscent of patterns observed in 
liquid flows including Gortler vortices in curved channel 
flow(ll], [[2) or longitudinal vortices in laminar/turbulent- 
boundary-layer transition [|l3| [l4|, |2lJ. However, the 
Gortler instability is driven by the inertial forces coming 
from the curvature, while streamwise vortices in bound- 
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ary layers seem to be linked to the turbulent nature of 
the flow Therefore, neither of these mechanisms 

can explain the structure we observed in our granular 
flow configuration. 

To understand the instability mechanism, we have 
studied how the volume fraction (the ratio between the 
volume occupied by the particles and the total volume, 
called simply density in the following) varies in the flow. 
The first result concerns the mean density measured by 
trapping some material during the flow and weighing the 
trapped material. The mean density is found to be never 
higher than 0.3 ± 0.05 in the whole range of parameters 
where the instability is observed. This result means that 
the instability takes place in a dilute flow regime. Sec- 
ondly, the density is not uniform across the bed when 
the flow becomes unstable but varies periodically with 
the surface deformation: crests are dilute and troughs 
are dense. Experimental evidence of this density vari- 
ation is presented in Fig. 4a. To obtain this picture, 
the transparent rough bed is lit from below and the light 
transmitted through the flow is collected from above by 
a CCD camera. The light transmitted by the thick part 
(the crest) is greater than the light transmitted by the 
thin part (the trough) . This is possible only if the crests 
are more dilute than the troughs. 

A complete picture of the flow pattern in a cross sec- 
tion can then be sketched in order to present the surface 
deformation, the grain motion and the density (Fig. 4b). 
One observes that the dense (thus heavy) part of the flow 
is going down while the dilute (thus light) part is going 
up, suggesting that density plays an important role in the 
instability. We thus propose the following mechanism to 
explain the formation of the granular longitudinal vor- 
tices. When the instability appears, the flow is rapid 
and dilute: its dynamics are controlled by the particle- 
particle and particle-boundary collisions. In this regime, 
the granular material can be seen as a dissipative dense 
gas, and a granular temperature can be defined related 
to the fluctuating motion of the grains |lq| . In our flow 
down a rough inclined plane, the source of the fluctuating 
motion is the bed roughness. Hence, as the flow acceler- 
ates from the outlet of the reservoir, particles close to the 
plane become more and more agitated due to collisions 
with the rough bed. The bottom granular temperature 
then increases along the slope. Consequently, the density 
at the bottom decreases (as in a molecular gas, density 
in a granular gas decreases when increasing the temper- 
ature) and eventually becomes smaller close to the plane 
than above. The flow is then mechanically unstable un- 
der gravity because the heavy material is above the light 
one yielding convective longitudinal rolls. There is a close 
analogy between this situation and the case of a liquid 
heated from below (Rayleigh-Benard instability). When 
a liquid film is flowing down an inclined hot plate, one 
observes longitudinal convective rolls j2^]. However, in 
this case the heating is provided by the imposed temper- 
ature at the bottom plate, whereas in the granular flow, 
the heating is produced by the flow itself. The rough sur- 



face induces a strong shear at the bottom which agitates 
the material from below. 
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FIG. 4: a, Picture obtained by lighting the flow from below 
to obtain information on the local volume fraction. The plots 
give the corresponding free surface transverse deformation, b, 
Sketch of the flow in a cross-section. 



The mechanism we propose is based on the density pro- 
file inversion. In a granular dissipative gas, the density 
profile actually results from a complex balance between 
gravity, collisions and dissipation and its prediction is not 
straightforward. Density profiles with higher density at 
the free surface than below have been indeed observed 
in numerical simulations of rapid flows using discrete el- 
ement methods , 23 . However, no instability was ob- 
served because the simulations were two dimensional. In 



order to investigate the relevance of the proposed insta- 
bility mechanism, we have performed a linear stability 
analysis of steady uniform flows in the framework of the 
kinetic theory of rapid granular flows §1| |2~i| , |2q , |(| $7\ . 
We have found that inverted density profiles are un- 
stable to transverse perturbations and yield to longitu- 
dinal rolls with transverse variations of velocities and 
density in qualitative agreement with the experimental 
observations [p8|. 

The new instability presented in this paper leads to a 
spontaneous pattern formation in a granular flow. An 
important issue is whether the longitudinal vortices rep- 
resent the first step in a global scenario leading to more 
and more complex structures. This evolution towards 
disordered states is well known in fluid mechanics for the 
transition to turbulence, but could be strongly affected 
in granular flows by the dissipative nature of the grain 
interactions. 
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